# Load data
nsduh <- read.csv("../Data/Modified_2022_NSDUH.csv", stringsAsFactors = TRUE)

# Subset youth only
youth <- nsduh[nsduh$is_minor, ] # subset for all minors
rm(nsduh) # dump original df

# Regressions
paratt.padsmoker.assmod <- lm(risk.pad_cigs ~ youth.parattitude_padsmoker, data = youth)
paratt.dailydrink.assmod <- lm(risk.bingedrink_daily ~ youth.parattitude_dailydrink, data = youth)
diff.herointry.assmod <- lm(risk.heroin_try ~ predictor.difficulty_heroin, data = youth)
diff.cocainemonthlyuse.assmod <- lm(risk.cocaine_monthlyuse ~ predictor.difficulty_cocaine, data = youth)

# Export models
save(paratt.padsmoker.assmod,
     paratt.dailydrink.assmod,
     diff.herointry.assmod,
     diff.cocainemonthlyuse.assmod,
     file = "../Results/Analysis_2.rdata")